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The idea of building a brain-inspired cognitive system has been around for several decades. Recently, 
electric-double-layer transistors gated by ion conducting electrolytes were reported as the promising 
candidates for synaptic electronics and neuromorphic system. In this letter, indium-zinc-oxide transis- 
tors gated by proton conducting methylcellulose electrolyte films were experimentally demonstrated 
with synaptic plasticity including paired-pulse facilitation and spatiotemporal-correlated dynamic 
logic. More importantly, a model based on proton-related electric-double-layer modulation and 
stretched-exponential decay function was proposed, and the theoretical results are in good agreement 
with the experimentally measured synaptic behaviors. © 2016 AIP Publishing LLC. 


[http://dx.doi.org/10.1063/1.4941080] 


Conventional computers based on CMOS technology and 
von Neumann architecture excel in logic and accurate scien- 
tific calculations, but it cannot compete with biological brains 
where the problems involve interaction with the real world and 
where both the inputs and outputs are often imprecisely speci- 
fied.' Inspired by the energy-efficient and cognitive ability of a 
biological brain, neuromorphic computing has emerged as an 
attractive computation paradigm that complements the von 
Neumann architecture.** In order to realize neuromorphic 
computing, software- and hardware-based approaches have 
been proposed. However, the software-based approaches are 
faced with significant challenges in energy dissipation since 
the algorithm is essentially run by conventional von Neumann 
computers with limited parallelism. The hardware-based 
approaches can potentially overcome such challenges by 
implementing massive parallelism at a physical level. Synaptic 
electronics have aroused wide attention in recent years, which 
aim at developing artificial synaptic devices to emulate the 
memory and computation performed by biological synapses. 
Up until now, a broad spectrum of synaptic devices has been 
developed with essential synaptic functions. !! Emerging 
memory structures have several characteristics such as 
historical-dependent analogue-like state, which are suitable for 
neuromorphic systems. For example, amorphous InGaZnO 
memristor was proposed to mimic an experimentally observed 
form of Hebbian learning rule known as spike-timing-depend- 
ent plasticity (STDP), which is mainly due to oxygen ion 
migration/diffusion induced conductance changes.'! 

Carriers in semiconductor channel layers can be strongly 
accumulated or depleted by electrolyte gating due to the 
extremely large electric-double-layer (EDL) capacitance formed 
at the electrolyte/semiconductor interface. !4 The EDL can be 
regarded as a nanoscale capacitor with a huge specific capaci- 
tance above 1.0 uF/cm?. Recently, three-terminal EDL transistors 
gated by ion conducting electrolytes were also proposed as the 
promising candidates for synaptic electronics.” "° In previous 
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reports, the gate pulse and channel current (or channel conduct- 
ance) were defined as the synaptic stimulus and synaptic weight, 
respectively. The synaptic weight of the transistors can thus be 
modified by tuning the pattern and timing of the stimuli, which 
enable the plasticity behaviors as a biological synapse. However, 
a physical insight into such synaptic devices is still lacking. 

Indium zinc oxide (IZO) film with high mobility 
(10-50 cm*/Vs) and large-area uniformity can be deposited by 
radio-frequency sputtering at a low temperature or even at a 
room temperature.”” Thin-film transistors with IZO channel 
layers are very good candidates for switching or driving tran- 
sistors in the next generation and/or flexible flat panel displays, 
sensors, and other new-concept applications, such as neuro- 
morphic systems.”!* In this letter, IZO-based EDL transistors 
gated by proton conducting methylcellulose electrolyte films 
were fabricated and used for synaptic function emulation. 
Synaptic plasticity emulations were also theoretically modeled 
based on the proton-related EDL modulation and stretched- 
exponential decay function, and the theoretical results were in 
good agreement with the experimental results. These findings 
could provide a useful guideline for designing artificial synap- 
ses for neuromorphic computing systems. 

Methylcellulose (MC) is a chemical compound derived 
from cellulose, which has good thermal, chemical, mechanical 
stabilities. It is synthetically produced by a substitution reac- 
tion that the hydroxyl residues (-OH functional groups) are 
replaced by methoxide (-OCH3 groups). These functional 
groups enable the MC film to absorb water molecules, which 
contributes to the high proton conductivity.” MC film with a 
good stability and a high proton conductivity is of great inter- 
est for EDL transistors’ fabrication and synaptic electronics 
application. First, MC powder (Sinopharm Chemical Reagent 
Co., Ltd.) was mixed into deionized water at 85°C for 5min 
and then cooled in air ambient to obtain a homogeneous solu- 
tion with a concentration of 1.0 wt. %. Then, the solution was 
spin-coated onto the ITO glass substrate and dried in air ambi- 
ent at 50°C for 2h. The thickness (D) of the MC film is esti- 
mated to be ~9.0 um. After that, 150-nm-thick IZO patterns 
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(1000 um x 150 um) were deposited on the MC film by mag- 
netron sputtering at room temperature with the use of a nickel 
shadow mask. As shown in Fig. 2(a), when the gap between 
the shadow mask and the substrate is about 50 um, a thin IZO 
channel layer with the thinnest part of ~30nm can be self- 
assembled between the IZO source/drain electrodes due to the 
reflection of IZO nanoparticles at the mask edge and dimen- 
sional extension by IZO nanoparticles with a low incident 
angle.'’** Proton conductivity and specific capacitance of the 
MC films were measured by a Solartron 1260A impedance an- 
alyzer. The electrical performances of the IZO synaptic tran- 
sistors were characterized by using a Keithley 4200 SCS 
semiconductor parameter characterization system. 

The proton conductivity of the MC film is determined 
from the Cole-Cole plot, as shown in Fig. 1(a). Proton conduc- 
tivity (a) is calculated by a relation o = D/(R — Ro)A, where 
D, R, Ro, and A are the thickness of MC film, the intersecting 
point on the real axis of Cole-Cole plots, the resistance of the 
electrodes, and the resistance of the electrode area, respec- 
tively. Here, A is 1.5 x 10-*cm’, and Ro is ~30Q. 
Therefore, the o value is estimated to be ~6.5 x 10°* S/cm. 
Figure 1(b) shows the frequency dependent specific capaci- 
tance of the proton conducting MC film measured with an 
in-plane test structure (inset in Fig. 1(b)). The specific capaci- 
tance increases with decreasing frequency and gradually 
reaches the maximum value (Cp) of ~3.1 uF /cm* at 1.0 Hz, 
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FIG. 1. (a) The proton conducting characteristics of the methylcellulose elec- 
trolyte film: Cole-cole plot. Inset: the molecular formula of methylcellulose. 
(b) The frequency dependent specific capacitance of the methylcellulose 
electrolyte film. Inset: the capacitance test structure with two in-plane IZO 
electrodes. 


which indicates the formation of EDL at the electrode/electro- 
lyte interface.’ 

Figure 2(a) shows the schematic image of IZO-based 
synaptic transistor gated by MC film on the glass substrate 
with two in-plane gates. The voltage applied on the in-plane 
gates can be capacitively coupled to the IZO channel.”* 
Figure 2(b) shows the transfer curve of the IZO-based synap- 
tic transistor with a constant Vps of 1.5 V. The gate voltage 
was swept from — 1.0 V to 0.8 V and then back. A clear anti- 
clockwise hysteresis window of ~0.25V was observed, 
which is likely due to the mobile protons within the MC 
electrolyte film.'* The subthreshold swing (SS), current on/ 
off ratio (Ipn/Iorr), and the threshold voltage (Vrp) are esti- 
mated to be ~85mV/decade, ~2.5 x 10°, and ~0.14 V, 
respectively. In addition, the field-effect electron mobility 
(ure) at the saturation region (Vps > Ves — Vrp) is esti- 
mated to be ~21.2 cm?/V.s by the following equation: 
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where W, L, and Cp, are the channel width, channel length, 
and specific capacitance of the EDL, respectively. 

Excitatory postsynaptic currents (EPSCs) or inhibitory 
PSCs (IPSCs) can be triggered by the potential spikes from 
the presynaptic neurons, which are collectively processed by 
the postsynaptic neurons to establish spatial and temporal 
correlated functions.” As previously reported, the gate pulse 
and channel current (or channel conductance) were defined 
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FIG. 2. (a) Schematic image of the fabrication of IZO-based synaptic tran- 
sistor gated by methylcellulose electrolyte film by a nickel shadow mask 
self-assembling method. (b) Transfer characteristics of the IZO-based synap- 
tic transistor operated at Vps = 1.5 V. The square root of Ips vs Vas is 
plotted as the blue line. 
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as the synaptic stimulus and synaptic weight, respectively. 
The synaptic weight of the devices can be modified by tuning 
the pattern and timing of the stimuli, which is similar to that 
in a biological synapse. In our experiment, presynaptic 
spikes were applied on the gate electrodes in order to emu- 
late the synaptic behaviors. The channel current measured at 
Vps =0.1 V was regarded as EPSC. Figure 3(a) shows the 
EPSC response triggered by a presynaptic spike (40 ms, 
0.1 V) applied on Gatel (G1). The inset of Fig. 3(a) shows 
the schematic diagram for EPSC measurement. As shown in 
Fig. 3(b), during a presynaptic spike, protons are electrically 
triggered and migrated to the IZO channel/methylcellulose 
electrolyte interface. The transient accumulation of protons 
then induces a sharp increase in the channel current by the 
strong EDL electrostatic coupling. When the spike ends, pro- 
tons will diffuse back to the equilibrium positions due to the 
concentration gradient.” The decay of proton concentration 
at the IZO/methylcellulose interface thus results in a gradual 
decay of the channel current. 

As mentioned above, the increment of the channel cur- 
rent is due to the accumulation of protons at the MC/IZO 
interface. We assumed that there is no interfacial electro- 
chemical reaction at a low gate voltage (Vas < 1.5 y" a 
Thus, the electrolyte film can be simply modeled as a 
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FIG. 3. (a) A typical EPSC response triggered by a presynaptic spike 
(40 ms, 0.1 V) applied on the in-plane gate electrode. The experimental data 
are plotted as open circles. The fitting curve is plotted as a solid line. Inset: 
the schematic illustration of the measurements of EPSC. (b) Schematic dia- 
gram illustrating the dynamic processes of an EPSC event. The EDL can be 
equivalent to a resistor-capacitor (RC) circuit with a resistor (Rp) and a 
capacitor (Cp) connected in parallel. 


resistor-capacitor (RC) circuit, as shown in Fig. 3(b). The 
voltage across the electrolyte film yields 


VpL(t) = Va - i — exp (- malt (2) 


where V4 and Tp, are the applied voltage and time constant 
of the RC circuit, respectively. During the charging process, 
the induced carriers in the IZO channel layer will increase 
the channel conductance. If the duration of a positive spike 
is T, the proton concentration (t < T) will increase with the 
charging voltage. According to Eq. (2), the proton concentra- 
tion (AN(t), < r) can be written by the following equation: 


t 
AN(t),<7¢ = CpoL(Ves — Vry) i — exp (- =)| | Z4WLòr, 


(3) 


where Z, q, and ôg are the ionic valence (Z= 1 for proton), 
elementary charge, and the depth of the ion-accumulation 
layer in electrolyte, respectively. When the spike ends, pro- 
tons will diffuse back to the equilibrium position gradually. 
Stretched-exponential decay has been widely applied for mod- 
eling relaxation process such as hydrogen diffusion in «-Si:H, 
oxygen vacancies’ diffusion in the WO,, etc.?™> In our case, 
the stretched-exponential function can also be used for model- 
ing the proton relaxation processes in the proton conducting 
MC film. Thus, the proton concentration after the presynaptic 
spike (t > T) can be determined by the following equation: 


AN. = AN(T) cr exp p (=) ] O 


where qt and f are the time constant of ion relaxation and the 
stretch index ranged between 0 and 1, respectively. It should 
be noted here that the total EPSC (Ips) can be determined by 
the relation: Ips = Jp + Alps, where Ip is the base channel 
current when no presynaptic spike is triggered (Vgs =0 V), 
and Alps is the changes in the channel current (EPSC) 
induced by a presynaptic spike. The changes in conductivity 
(Ao) of the semiconductor layer could be determined by the 


following relations:*** 
AJ AlpsbL 
Ao = — =-—§ (5) 
E VpsWos 
and 
Ao = Anqu, (6) 


where AJ, E, ôs, and An are the changes in current density, the 
electric-field intensity, the depth of the carrier-accumulation 
layer in the semiconductor, and the induced electron concentra- 
tion, respectively. Therefore, AIpg can be written as 


Anqu - VpsW ôs 


Alps = Z 


(7) 
Due to the huge EDL capacitance, the induced electron con- 
centration should be approximately equal to the induced proton 
concentration (An = AN). Thus, the EPSC amplitude (A(t)) can 
be deduced by substituting AN(t), < r and AN(t), > q into Eq. (7) 
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where « = Cròs ee ). k can be regarded as the EPSC 


amplitude triggered by a continuous presynaptic spike (T — 00). 
The measured EPSC is plotted as open circles in Fig. 3(a). 
The fitted result based on Eq. (8) is plotted as a solid line in 
Fig. 3(a). The fitting parameters of K, ToL, T, and f are esti- 
mated to be 766nA, 257 ms, 9.0 ms, and 0.74, respectively. 
The fitting curve approaches the experimental results very 
well, indicating the validity of our proposed model. 
Paired-pulse facilitation (PPF) is a short-term form of 
synaptic plasticity, which represents an amplitude increase 
in the second synaptic response when it follows shortly after 
the first one.°°°7 Figure 4(a) shows the EPSC response trig- 
gered by the two successive presynaptic spikes (0.1 V, 
10 ms) with a time interval (AT) of 40 ms. The second EPSC 
peak (A2) is ~156nA, which is larger than the first EPSC 
peak (A,) of ~148 nA. The inset in Fig. 4(a) schematically 
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FIG. 4. (a) EPSCs triggered by two successive presynaptic spikes (0.1 V, 
10 ms) with a time interval (AT) of 40 ms. The black open circles are the 
experimental data, and the blue solid line is the fitting curve. Inset: the 
schematic diagram for measuring PPF. (b) The facilitation ratios (A2/A,) 
are plotted as a function of AT. The black open circles are the experimental 
data, and the blue solid line is the fitting curve. 


illustrates the measurements for PPF emulation. The facili- 
tation ratios (A2/A,) are plotted as a function of AT, as 
shown in Fig. 4(b). The obtained PPF emulation results 
are in good agreement with that observed in biological 
systems.°°°* According to Eq. (8), the EPSC triggered 
by the first spike would reach its maximum value (A,) of 
Kk- [1 — exp(— =, )] at the end of the spike (t=T). After 
that, the protons will diffuse back to their equilibrium posi- 
tions yielding the stretched exponential function. During 
the second pulse, the residual protons triggered by the first 
spike will be added to the total proton concentration, which 
would induce an increment in EPSC as a result. Therefore, 
the amplitude of the second EPSC peak can be estimated by 
the following relation: 


=x [t-e0(-Z)] e [t-en(-2)] 
ales] > 


where / is a factor related to the effective residual proton 
concentration contributed to the second EPSC. Then, the 
facilitation ratio (A2/A,) between the two EPSC amplitudes 
can be determined by the following equation: 


1 
The facilitation ratio (A2/A,) is plotted as a function of AT 
in Fig. 4(b). The blue solid line is the fitted by Eq. (10) with 
A=2.2. We find that the fitting curve approaches the experi- 
mental results very well. 

Dynamic logics established by the spatiotemporal cor- 
related spikes are also demonstrated in the proposed synap- 
tic transistor with two in-plane gates (Gl and G2) as 
presynaptic terminals. Figure 5(a) shows the schematic dia- 
gram of such an emulation. Two presynaptic spikes P1 and 
P2 (0.1 V, 10 ms) with a time interval (AT) are applied on 
two presynaptic terminals Gl and G2, respectively. The 
zero time (t=0) is defined as the time when P2 ends. The 
EPSC amplitudes measured at t = 0 are plotted as a function 
of AT, as shown in Fig. 5(b), which decrease asymmetri- 
cally with increasing |AT|. The EPSC amplitudes decrease 
gradually from 520nA to 217nA when AT changes from 0 
to —150 ms, while no significant change in the EPSC ampli- 
tudes can be observed when AT > 0. 

When AT > 0, i.e., the spike P1 triggers after spike P2, 
the EPSC amplitude measured at t=0 is equal to the EPSC 
individually triggered by P2. When AT <0, i.e., the spike P1 
triggers before spike P2, the EPSC amplitude measured at 
t=0 is the summed current of EPSC1 triggered by P1 and 
EPSC2 triggered by P2. Therefore, the summed current can 


B 
pian) == 142. ep|-(#2) | (10) 


be written as /’A, - exp[—(=42)4] + A>, where 2! is a factor 
related to the effective residual proton concentration (trig- 
gered by P1) contributed to the summed current, and A; = K1 
[1 — exp(— 1 )] and Ay = kp - [1 — exp(—)] are the EPSC 
amplitudes individually triggered by P1 and P2, respectively. 
In conclusion, the EPSC amplitude at t= 0 can be determined 
by the following equation: 
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FIG. 5. (a) Schematic diagram of the measurement for dynamic logic estab- 
lished by spatiotemporal correlated spikes. (b) The EPSC amplitude measured 
at t=0 is plotted as a function of time interval (AT). The experimental results 
are plotted as black open circles, and the blue solid line is the fitting curve. 


T B 
A'A: + exp (=) +A, AT<O 
T 
Ad AT > 0. 
(11) 


The EPSC amplitudes measured at t=O are plotted as a 
function of time interval (AT) in Fig. 5(b). Furthermore, the 
measured EPSC amplitudes are fitted by Eq. (11) and plotted 
as blue solid line in Fig. 5(b). 4’, Aj, and A; are estimated to 
0.92, 332nA, and 217nA, respectively. Again, we find that 
the fitting curve approaches the experimental results very 
well, indicating the validity of the proposed model. 

In summary, synaptic plasticity including PPF and 
spatiotemporal-correlated dynamic logic were experimentally 
demonstrated in an organic/inorganic hybrid artificial synapse 
based on an IZO-based EDL transistor gated by methylcellu- 
lose electrolyte film. A model based on proton-related electric- 
double-layer modulation and stretched-exponential decay 
function was proposed, and the theoretical simulation results 
are consistent with the measured synaptic behaviors. Our 
results may provide an interesting guideline for artificial syn- 
apses and brain-inspired neuromorphic systems. 
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